Dephasing of a non-relativistic quantum particle 
due to a conformally fluctuating spacetime 



Paolo M. Bonifacio 1 , Charles H.-T. Wang 1 2 , J. Tito 
Mendonga 2 3 , Robert Bingham 2,4 

X SUPA Department of Physics, University of Aberdeen, King's College, 
Aberdeen AB24 3UE, UK 

2 Rutherford Appleton Laboratory, STFC, Chilton, Didcot, Oxfordshire OX11 
OQX, UK 

3 CFP and CFIF, Institute Superior Tccnico, 1049-001 Lisboa, Portugal 
4 SUPA Department of Physics, University of Strathclyde, Glasgow G4 ONG, UK 

E-mail: p.bonif acioOabdn. ac .uk, cwangOabdn.ac.uk, titomendOist.utl.pt, 
r . binghamOr 1 .ac.uk 

Abstract. We investigate the dephasing suffered by a nonrelativistic quantum 
particle within a conformally fluctuating spacetime geometry. Starting from 
a minimally coupled massive Klein-Gordon field, the low velocity limit yields 
an effective Schrodinger equation where the wave function couples to gravity 
through an effective nonlinear potential induced by the conformal fluctuations. 
The quantum evolution is studied through a Dyson expansion scheme up to 
second order. We show that only the nonlinear part of the potential can induce 
dephasing. This happens through an exponential decay of the off diagonal terms 
of the particle density matrix. The bath of conformal radiation is modeled in 
3-dimensions and its statistical properties arc described in general in terms of a 
power spectral density. The case of a Lorcntz invariant spectral density, allowing 
to model vacuum fluctuations at a low energy domain, is investigated and a general 
formula describing the loss of coherence derived. This depends quadratically on 
the particle mass and on the inverse cube of a typical particle dependent cutoff 
scale. Finally, the possibilities for experimental verification are discussed. It 
is shown that current interferometry experiments cannot detect such an effect. 
However this conclusion may improve by using high mass entangled quantum 
states. 



PACS numbers: 03.65.-w, 03.65. Yz, 04.20.Cv, 05.40.-a 

1. Introduction 

It is generally agreed that the underlying quantum nature of gravity implies that the 
spacetime structure close to the Planck scale departs from that predicted by General 
Relativity. Unfortunately the Quantum Gravity domain is still beyond modern particle 
accelerators such as LHC. Nonetheless, finding experimental ways to test the quantum 
structure of spacetime would be highly beneficial to the theoretical developments of our 
fundamental theories of nature. In this respect it is has been suggested that quantum 
gravity could induce decoherence on a quantum particle through its underlying Planck 
scale spacetime fluctuations [1, 2, 3, 4, 5, 6]. 

As the sensitivity and performance of matter wave interferometers is increasing 
[7, 8, 9, 10, 11], it is important to assess the theoretical possibility of a future 
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experimental detection of intrinsic, spacctime induced decoherence. The closely 
related dephasing effect due to a random bath of classical GWs (e.g. of astrophysical 
origin) has been extensively studied e.g. in [12]. The problem of the decoherence 
induced by spacctime fluctuations is difficult to study as a quantum gravity theory 
is still missing. Notwithstanding promising progress, mainly in loop quantum 
gravity and superstring theory [13], a coherent and established quantum gravity 
theoretical framework is still missing. Thus any theoretical attempt for a prediction 
of the decoherence induced by spacetime fluctuations must exploit some semiclassical 
framework. Such approaches typically represent the spacetime metric close to the 
Planck scale by means of fluctuating functions. These are usually supposed to mimic 
the vacuum quantum property of spacetime down to some cutoff scale i = \Lp, where 
Lp is the Planck scale. The adimensional parameter A marks the benchmark between 
the fully quantum regime and the scale where the classical properties of spacetime 
start to emerge [1] . The fact that classical fluctuating fields can be used to reproduce 
various genuine quantum effects is well known, e.g. from the work of B oyer [14, 15] in 
the case of the EM field or Frederick [16] in the case of spacetime fluctuations. This 
is often exploited in the literature in relation to problems involving the microscopic 
behavior of the spacetime metric; e.g. a stochastic metric was employed in [17, 18] 
to study the problem of gravitational collapse and big bang singularities, while in 
[19] spacetime metric fluctuations were introduced and their ability to induce a WEP 
violation studied. 

A pioneering analysis of the problem of spacetime induced decoherence has been 
proposed by Power & Percival (PP in the following) [1] in the case of a conformally 
modulated Minkowski spacetime with conformal fluctuations traveling along 1 space 
dimension. This was improved by Wang et al. [3], who extended upon PP work 
attempting to include the effect of GWs. Conformal fluctuations are interesting as 
they are mathematically easy to treat and offer a convenient way to build 'toy' models 
to assess some of the problem's features. They have an important role in theoretical 
physics [20] and are sometimes invoked in the literature also in relation to universal 
scalar fields [21, 18] that can arise naturally in some modified theories of gravity such 
as scalar-tensor theories [22, 23]. 

Within a semiclassical approach that 'replaces' the true quantum environment by 
classical fluctuating fields we should properly speak of dephasing of the quantum 
particle rather than decoherence. In a remarkable paper [24] about quantum 
interference in the presence of an environment, Stern and co-authors showed that 
a fully quantum approach that studies decoherence by 'tracing away' the environment 
degrees of freedom in the quantum system made up by system + environment, and 
that in which the dephasing of the quantum particle is due to a stochastic background 
field give equivalent results. 

In this paper we consider a conformally modulated 4-dimcnsional spacetime 
metric of the form g a b = (1 + A) 2 n a i,. Such a metric has been considered by PP 
[1] , where the dephasing problem was studied in the simple idealized case of a particle 
propagating in 1-dimension. By imposing Einstein's equation on the metric g a b, PP 
deduced a wave equation for A. Their procedure to derive an effective newtonian 
potential interacting with the quantum probe started from the geodesic equation of 
a test particle. Even though this didn't take properly into account the nonlincarity 
in the conformal factor (1 + A) 2 , they found correctly that the change in the density 
matrix is given by Sp oc M 2 TAqT*, where M is the probing particle mass, T the flight 
time, Aq the amplitude of the conformal fluctuations and their correlation time. 
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This formula was used to set limits upon A. However in doing this they did not treat 
the statistical properties of the fluctuations properly and this resulted in the wrong 
estimate A oc (M 2 T '/ 'dp) 1 / 1 ', as already noted by Wang et al in [3]. 

In their work, Wang and co-authors attempt to include GWs into the analysis 
by considering a metric of the kind g a b = (1 + A) 2 ^^. This was done by 
exploiting the results in [25, 26] where a canonical gcomctrodynamics approach 
employing a conformal spacial 3-metric was studied. By exploiting an energy density 
balancing mechanism between the conformal and GWs parts of the total gravitational 
Hamiltonian the statistical properties of the conformal fluctuations where fixed. This 
corresponded to assume that each 'quantum' of the conformal field possessed a zero 
point energy —hw. Though an improvement over PP work, this approach is still 1- 
dimensional and too crude to make predictions. Moreover, as it shall be discussed 
extensively in a future report [27], the issue of energy balance between conformal 
fluctuations and GWs is a delicate one, and likely not to occur within the standard 
GR framework. 

In the present work we provide a coherent 3-dimensional treatment of the problem 
of a slow massive test particle coupled to a conformally fluctuating spacetime. The 
conformal field A is assumed to satisfy a simple wave equation. This will allow a direct 
comparison with PP result. We also notice that such a framework is expected to arise 
naturally within a scalar-tensor theory of gravity. This issue will be discussed in a 
future report [28]. 

The work is organized as follows: in section 2 the correct non-relativistic limit 
of a minimally coupled Klein-Gordon field is deduced and an effective newtonian 
potential depending nonlinearly on A is identified in the resulting effective Schrodinger 
equation. In section 3 we set the general formalism to study the average quantum 
evolution through a Dyson expansion scheme for the particle density matrix p. In 
section 4, general results derived in Appendix A are used to model the statistical and 
correlation properties of the fluctuations through a general, unspecified, power spectral 
density. In section 5 and 6 wc compute the average quantum evolution and derive a 
general expression for the evolved density matrix. We show in general that only a 
nonlinear potential can induce dephasing. The resulting dephasing formula implies an 
exponential decay of the density matrix off-diagonal elements and is shown to hold in 
general and independently of the specific spectral properties of the fluctuations. All 
we assume is that these obey a simple wave equation and that they are a zero mean 
random process. The overall dephasing predicted within the present 3-dimensional 
model -equation (28)- is seen to be about two orders of magnitudes larger than in 
the 1-dimcnsional case as derived by PP. Next we consider in section 7 the problem 
of vacuum fluctuations. To this end a power spectrum S(uj) ~ l/u> is introduced 
and we derive an explicit formula for the rate of change of the density matrix. This 
result improves over both Percival's and Wang's work in that its key ingredients are 
general enough to be potentially suited for a variety of physical situations. Finally 
the discussion in section 8 addresses the question of whether the dephasing due to 
conformal vacuum spacetime fluctuations could be detected. In other words, whether 
the proposed theory can be falsified or not. A possibility would be through matter 
wave interferometry employing large molecules. We consider this issue in the final part 
of this paper by estimating the probing particle resolution scale, setting its ability to 
be affected by the fluctuations. The resulting formula for the dephasing rate indicates 
that the level of the effect is still likely to be beyond experimental capability, even 
for large molecules (e.g. fullerenes) [10]. A measurable effect could possibly result for 
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larger masses, e.g. if entangles quantum states were employed [29]. 

2. Low velocity limit and effective Schrodinger equation 

The problem we wish to solve is clearly defined: we consider a scalar field A inducing 
conformal fluctuations on an otherwise flat spacetime geometry according to 



where r\ a b = diag(— 1, 1, 1, 1, ) is the Minkowski tensor. We will refer to A as to 
the conformal field and this will be assumed to satisfy the wave equation d c d c A = 
0. Solving this equation with random boundary conditions results in a randomly 
fluctuating field propagating in 3-dimcnsional space. We assume this to be a small first 
order quantity, i.e. \A\ — 0(e <C 1). Equation (1) expresses the spacetime metric in 
the laboratory frame. We also suppose that the typical wavelengths of A are effectively 
cut off at a scale set by I := ALp, where Lp = (hG/c 3 ) 1 ^ 2 ~ fO -35 m is the Planck 
length. The adimcnsional parameter A represents a structural property of spacetime 
marking the quantum-classical transition: below £ a full quantum treatment of gravity 
would be needed so that, by definition, I represents the scale at which a semiclassical 
approach that treats quantum effects by means of classical randomly fluctuating fields 
is supposed to be a valid approximation. The value for A is model dependent but it is 
generally agreed that A > fO 2 [3], so that I is expected to be extremely small from a 
macroscopic point of view. This motivates the assumption that classical macroscopic 
bodies, including the objects making up the laboratory frame and also the observers, 
arc unaffected by the fluctuations in A. This corresponds to the idea that a physical 
object is characterized by some typical resolution scale Lr that sets its ability to 'feel' 
the fluctuations: if Lr, 3> I these average out and do not affect the body, that simply 
follows the geodesic of the flat background metric. On the other hand a microscopic 
particle can represent a successful probe of the conformal fluctuations if its resolution 
scale is small enough. 

We are interested in the change of the phase induced on the wave function of 
a quantum particle by the fluctuating gravitational field. Various approaches to the 
problem of how spacetime curvature affects the propagation of a quantum wave exist 
in the literature; e.g. for a stationary, weak field and a non relativistic particle a 
Schrodinger -like equation can be recovered [30]. The more interesting case of time 
varying gravitational fields can be treated e.g. by eikonal methods that are usually 
restricted to weak fields with g a b = n a b + h a b and \h a b\ <C 1 [3f, 32]. Other approaches, 
e.g. in [f9, 6], are based on the scheme developed by Kiefer [33] for the nonrelativistic 
reduction of a Klein-Gordon field which is minimally coupled to a linearly perturbed 
metric. The approach of PP in [1] and of Wang at el. in [3] was to derive the 
geodesic equation in the weak field limit. Their treatments were however employing, 
incorrectly, the usual newtonian limit scheme which is valid only for weak, linear and 
static perturbations [34]. Conceptually the wave approach is more satisfying than 
that based on the geodesic equation because, e.g., the coupling between gravity and a 
scalar field is well understood and in the appropriate non-relativistic weak field limit 
an effective Schrodinger equation emerges. This will be our approach below. 



We describe the quantum particle of mass M by means of a minimally coupled 
Klein-Gordon (KG) field <f>: 



9ab = (1 + A) 2 T]ab, 



(1) 



5 ab V a V h 
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where V a is the covariant derivative of the physical metric g a \,. Using g ab = f2 2 77 Q {, 
this equation can easily be made explicit [34] and reads: 

1 FP- \ Q 2 M 2 r 2 

~ + V 2 ) ^ - 2d a Qn 0)0"*, (2) 

i.e. the wave equation for a massive scalar field plus a perturbation due to A describing 
the coupling to the conformally fluctuating spacetime. We need to take the appropriate 
non-relativistic limit in order to deduce an effective Schrodinger equation. Before 
doing this we remark that, had we consider the alternative meaningful scenario of a 
conformally coupled scalar field, then the equation g ab \7 a \7 b (f)— Rcf)/6 — M 2 c 2 (j)/h 2 = 
would read explicitly 

1 8 2 \ fl 2 M 2 r 2 

+ y2 J ^ = ~1T 2 — ~ ^(^n)^ - W Weft. 

Since it is Q,~ 1 d c d c £l = (1 — A)d c d c A + 0(e 3 ) we see that if A is assumed to satisfy 
the wave equation then the curvature term has no effect: in this case the minimally 
and conformally coupled KG equations are equivalent up to second order in A. We 
also note that, by introducing the auxiliary field $ := f2</>, equation (2) turns out to 
be equivalent to: 

1 d 2 „ 2 \ Q 2 M 2 c 2 

hV $= 

c 2 dt 2 J h 2 

In principle, if a solution for $ were known, then the physical scalar field representing 
the particle would follow formally as <\> — = (1 — A + A 2 )$, up to second order in 

A. However in studying the dephasing problem we will only find an averaged solution 
for the average density matrix representing the quantum particle. Therefore, even 
if a solution in this sense is know in relation to it would not be obvious how to 
obtain the corresponding averaged density matrix related to <p, which is what we are 
interested in. 

In view of the above considerations we work directly with equation (2) and now 
proceed in deriving its suitable non-relativistic limit. We will make two assumptions: 

(i) the particle is slow i.e., if p = Mv is its momentum in the laboratory, we have: 

- < l; 

c 

(ii) the effect of the conformal fluctuations is small, i.e. the induced change in 
momentum Sp = MSv is small compared to Mv: 

Sv 

— < 1. 

v 

In view of these assumptions we can write: 

<j) = f/>exp (-iMc 2 t/h), 

where the field ip is close to be a plane wave of momentum p. As a consequence we 
have: 
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d 



dt 2M 



h 2 2 Mc 2 
V 2 — 



* A2 



T 2 



Mc 2 ip 



5* 



% ln(l + A) x cxp (iMc 2 t/ti) 



(3) 



T 3 



Leaving the term T4 aside for the moment, the orders of the three underlined 
terms must be carefully assessed. We have: 

4 



Mv 2 , 



To 



[ - 



e 2 Mc 2 



where an average has been inserted since T3 is fluctuating. It follows that 
To 



5f~l-l «i 

Thus, in the non-relativistic limit, T2 is negligible in comparison to T\. This is the 
case in a typical interferometry experiment where it can be v ~ 10 2 m s _1 [11], so 
that (v/c) 2 - 1CT 12 . Next we have: 



(Ts) 
Ti 



V v 



The request that the conformal fluctuations have a small effect thus gives the condition 



Ti 



< 1 <^ 



(A 2 ) « V - 



(4) 



That this condition is effectively satisfied can be checked a posteriori after the model 
is complete. It depends on the statistical properties of the conformal field and the 
particle ability to probe them. This will be related to a particle resolution scale. At 
the end of the discussion in section 8.3 we will show that (4) is satisfied if, e.g., the 
particle resolution scale is given by its Compton length. 

Under these conditions the non-relativistic limit of equation (3) yields: 



2M 



V 2 i/j 



A 



A 2 



Mc 2 tfj + T 4 = ih 



dt' 



(5) 



where 



Ti := -—d a <f>d a ln(l + A) X exp (iMc 2 t/ti). 



In order to assess the correction due to this term we split it into two contributions by 
writing separately the time and space derivatives. Using the fact that 

d<j> ^ iMc 2 
it is easy to see that: 



-■0exp (-iMc 2 t/h) 



-ih (A - AAj i' - ihv (A x - AA iX ) ip, 



(6) 
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where A := dA/dt, A x := dA/dx and where we assumed that the particle velocity 
in the laboratory is along the x axis. In Appendix B we show that if A is (i) a 
stochastic isotropic perturbation and (ii) effectively fast varying over a typical length 
\a = nh/(Mc) related to the particle resolution scale, then T4 reduces to: 

T 4 = (-A + A 2 )—i;. (7) 

K 

Here k ~ 1 is dimensionless and its precise value is unimportant. The important point 
is that T4 yields a positive extra nonlinear term in A that adds up to what we already 
have in (5). Finally we get the effective Schrddinger equation 

2M v v dt' 
where the nonlinear fluctuating potential V is defined by 

V := (dA + C 2 A 2 ) Mc 2 . (8) 
The values of the constants C\ and Ci depend on k and C2 is always strictly positive. 
For k = 1 it would be C\ = and C2 = 3/2. For generality we will leave them 
unspecified in the following treatment and consider k as a constant of order one. 



3. Average quantum evolution 

3.1. Dyson expansion for short evolution time 

We now have a rather well defined problem: that of the dynamics of a non 
relativistic quantum particle under the influence of the nonlinear potential (8). The 
Schrodingcr equation describing the dynamics of a free particle is suitable to describe 
the interference patterns that could result e.g. in an intcrferometry experiment 
employing cold molecular beams. When the particle in the beam propagates through 
an environment, we are dealing with an open quantum system. This in general suffers 
decoherence, resulting in a loss of visibility in the fringes pattern [9, 11]. This is a well 
defined macroscopic quantity. In the present semiclassical treatment the environment 
due to spacetime fluctuations is represented, down to the semiclassical scale £, by a 
sea of random radiation encoded in A and resulting in the fluctuating potential V. 
An estimate of the overall dephasing can be obtained by considering the statistical 
averaged dynamics of a single quantum particle interacting with V. In practise we will 
need (i) to solve for the dynamics of a single particle of mass M and (ii) calculate the 
averaged wavefunction by averaging over the fluctuations. The outcome of (i) would 
be some sort of 'fluctuating' wavefunction carrying, beyond the information related 
to the innate quantum behavior of the system, that related to the fluctuations in the 
potential. The outcome of (ii) is to yield a general statistical result describing what 
would be obtained in an experiment where many identical particles propagate through 
the same fluctuating potential. 

We thus consider the Hamiltonian operator H(t) = H° + H l {i) : where H° is the 
kinetic part while 

H\t) = J d 3 xV(x,t)\x)(x.\, 

is the perturbation due to the fluctuating potential energy. Here |x)(x| is the 
projection operator on the space spanned by the position operator eigenstate |x). 
Indicating the state vector at time t with ijj t , the related Schrodingcr equation reads 
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Using the density matrix formalism, the general solution can be expressed through a 
Dyson series as [36] 

p T = po + Ki(T) Po + p Kl(T) + K 2 (T)p + K^pokliT) + P J4(T) + ..., 
where po is the initial density matrix and the propagators K\ (T) and K 2 (T) are given 

by 

^i(T):=~^ H(t')dt', 

K 2 (T) := dt' J* dt"H(t')H{t"). 

In truncating the series to second order we assume that the system evolves for a time 
T such that T « T*, where T* is defined as the typical time scale required to have a 
significant change in the density matrix p. 

The effect of the environment upon a large collection of identically prepared 
systems is found by taking the average over the fluctuating potential as explained 
above. Formally and up to second order we have 

(pr) = (po + K x {T) Po + p Q K\{T) + k 2 (T)p Q + K^pokliT) + Po kl(T)) . 

The average density matrix ( P t) will describe the average evolution of the system 
including the effect of dephasing. 

It is straightforward to show that, up to second order in the Dyson's expansion, 
the kinetic and potential parts of the hamiltonian give independent, additive 
contributions to the average evolution of the density matrix, i.e. (pr) = [pt]o + ([pt]i) , 
where 

[prh ■■= Po + [£i(T)]opo + Po[kx(T)]l 

+ [k 2 (T)] po + [£x(T)]oPo[£i(T)]J + Po[k 2 (T)]l 

([Hi) := (po + [ki(T)) lPo + p Q [ki(T)]\ 

+ [k 2 (T)] lP0 + [£i(T)]ipo[£iCT)]J + P o[k 2 (T)}\) . (9) 

Here the kinetic propagators [iCi(T)] and [K 2 (T)] depend solely on H°, while the 
potential propagators [_K"i(T)]i and [k 2 (T))i depend only on In the next 

section we estimate the dephasing by calculating the term ([/3t]i) alone. 

4. The conformal field and its correlation properties 

We now set the statistical properties of the conformal field A. This is assumed 
to represent a real, stochastic process having a zero mean. We further assume it 
to be isotropic. In Appendix A we review a series of important results concerning 
stochastic processes, in particular in relation to real stochastic signals satisfying the 
wave equation. The main quantity characterizing the process is the power spectral 
density S(cu). In the case of an isotropic bath of random radiation, field averages such 
as (A 2 ), (\VA\ ) and ((d t A) 2 ) can be found in terms of S(u>), e.g. 

(A(x,tf) = ^Jd 3 kS(k), 
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where kc = u>. In Appendix A we show how the conformal field can be resolved into 
components traveling along all possible space directions according to 



A(x,t) = J dkA^k-x/c-t), 



where dk indicates the elementary solid angle. The capacity of the fluctuations to 
maintain correlation is encoded in the autocorrelation function C(r). In the same 
appendix we prove a generalization of the usual Wiener-Khintchinc theorem, valid for 
the case of a spacetime dependent process satisfying the wave equation, and linking 
the autocorrelation function to the Fourier transform of the power spectral density 
according to: 

C(t) := t-^t / duJLu 2 S(Lu)cos(LUT). (10) 

(2TTCy J 

This allows to prove that wave components traveling along independent space 
directions are uncorrelated, i.e. 

(A^t)A^(t + r)) = 6(k,-k')C(r). (11) 

The field mean squared amplitude is related to the correlation function according to 
(A 2 ) = 47rCo, as derived in Appendix A 

Isotropy implies that all directional components have the same amplitude Aq. 
This is found introducing the normalized correlation function R(t) through 

so that i?(0) = 1. Equation (11) can now be re-written as (A^t) A^,(t + r)) = 
£(k, k') Co R(t) so that, introducing the normalized directional components, ft(t) ■= 
A^it) / \/Co we have (f^(t) f^, (t + r)) = <5(k, k') i?(r). We now define the constant 
Aq :— ^/Cq which is connected to the squared amplitude per solid angle according to 
Aq = Co = (A 2 ) /4tt. The directional components are given by A^(t) = Aq f^(t) and 
the general conformal field can finally be expressed as an elementary superposition of 
the kind 

A(x, t)=A J die / fc (k ■ x/c - t). (12) 

4--1- Summary of the correlation properties of the conformal fluctuations 

The main statistical properties of the directional stochastic waves f^ are summarized 

by 

=0, (13) 
(ft(t)k,(t'))=5(k,k>)R(t-t'), (14) 

i.e. each component has zero mean and fluctuations traveling along different space 
directions are perfectly uncorrelated. These two properties imply that odd products 
of directional components have also a zero mean, i.e. 

(^i)^)/^)) =0. (15) 
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In the following dephasing calculation we will need to evaluate means involving 
products of four directional components. To this purpose we need to introduce the 
second order correlation function R"(t — t') according to 

([k(t)] 2 lk'(t')} 2 ) = 1 + S & - f ) - !]• (16) 

This definition is compatible with the fact that the mean is one when components 
traveling in different direction are involved, i.e. ([k{t)] 2 [fy(t')] 2 ) — 1 if k ^ k'. 



5. Dephasing calculation outline 



To calculate the dephasing suffered by the probing particle we must evaluate the 
average of all the individual terms in equation (9). The relevant propagators are 

rT 



ft Jo 



dt'H 1 ^'), 



[Kin ■■- 



[&2(T)]i ~ ~ ¥ I dt I dt'H'm^t'). 
The interaction Hamiltonian is given by 

H\t)= J d 3 xV(x,t)\x)(x\, 
where the potential energy is 

V(x, t) = dMc 2 A ( dk/ t (t - x • k/c) + C 2 Mc 2 A 2 



dk/ t (*-x-k/c) 



5.1. First order terms of the Dyson expansion 



We evaluate the two first order terms in the Dyson expansion. For a more compact 
notation, we do not show the argument of the directional components f^. The 
contribution of the linear part of the potential C\Mc 2 A vanishes trivially since 
(A) = 0. The quadratic part gives: 

2\ 



K l (T)p ) = - 



iC 2 Mc A A" 







j dt Jd 3 x\x)(x\ / 









Po- 



Using A^{t) = yCo4(i) and {A 2 ) = 4nCo = ^A\ it is seen that the average yields 
47T. Since §d?x |x)(x| = I and integrating over T we find 

A-KC 2 iMc 2 A 2 Q T 



ifi(T)A) 



- Pa- 



in) 



The calculation of the other first order term proceeds in the same way. Since 
K\{T) = —Kx, it yields the same result as in (17) but with the opposite sign (more 
in general, all the odd terms in the Dyson expansion have an i factor and also yield a 
vanishing contribution). We thus see that at first order in the Dyson expansion there 
is no net dephasing and (Ki(T)po + pqK\(T)) = 0. 
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5.2. Second order terms of the Dyson expansion 

The second order calculation is more complicated. A fundamental point is that the 
linear part of the potential does again give a vanishing contribution. Dephasing will 
be shown to come as a purely nonlinear effect due to the nonlinear potential term 
~ A 2 . 

5.2.1. (N on) -contribution of the linear part of the potential To have an idea of how 
things work we consider e.g. the average of the term K-iPo- This has the following 
structure: 

(K2Pv)~ JdtJ dt 1 J d 3 y|y)(y| J d 3 y'\y')(y'\ (V(y,t)V(y',t')) p . 

The interesting part is the average (V(y, i)V(y', t')). This is: 

(V(y,t)V(y',t')) ~ (A(y,t)A(y',t')) + (A(y,t)A 2 (y',t')) + (A 2 (y,t)A 2 (y',t')). 

The first two term are due to the linear part of the potential. The second of them 
vanishes in virtue of property (15). This is seen using the directional decomposition 
(12) and writing: 

(A(y,t)A 2 (y',t')) = A* J dk t J dk 2 J dk 3 (/ tl 4 2 4 3 ) = 0. 

The first term results in the contribution: 
(A(y,t)A(y',t')}^ 

[dtfdt'J d 3 y|y)(y| J d?y'\y'){y'\ J dk|dk' (/ fc (t - y • k)f p (t> - y' ■ k!)) Po- 

For convenience of notation we set c = 1 in the arguments of the directional functions 
fy.. Using equation (14) the average yields the 2-point correlation function according 
to S(k, U )R(t — t' + y' • k' — y • k). Integrating with respect to k yields: 

(A(y,t)A(y',t'))^ 

" dt jdt'j d 3 y|y)(y| J d 3 y'|y')(y'l J dkR{t - f + k ■ (y' - y)]p . 

The corresponding matrix element is found by inserting (x| and |x') respectively on 
the left and on the right. Using (x|y) = S(x — y) and exploiting the properties of the 
delta function we find 

(A(y, t)A{y',t')) => K x f dt [ dt' [ dkR(t - t'), (18) 

Jo Jo J 

where K is a constant given by 

CfA 2 M 2 c 4 p xx ,(0) 
K= h 2 ' 

and where p xx / (0) := (x|po|x'). The similar terms coming from (^paK^) will contribute 
in the same way as in (18), thus yielding an extra factor 2. Finally, through a similar 
calculation it is found that the terms ~ A(y,t)A(y',t') coming from (^KiPqK\ 
contribute according to: 

(A(y,t)A(y ! ,t')) -K x f dt f dt'[dkR[t - t' + k ■ (x' - x)]. 

Jo Jo J 
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Bringing all together, the overall contribution deriving from the linear part C\Mc 2 A 
of the effective potential is found to be proportional to the expression: 

I -=J o dt^2 J^dt'R{t-t')- J^dt'R[t-t' + k-(x' . (19) 

In Appendix C we prove that this vanishes provided R(t) is an even function and 
provided that the drift time T is much larger than the time needed by the fluctuations 
to propagate through the distance |Ax|, i.e. if T 3> k- Ax, where c = 1. This condition 
is certainly satisfied in a typical intcrfcromctry experiment where the drift time T 
can be of the order of ~ 1 ms and cT is indeed much larger than the typical space 
separations |Ax| relevant to quantify the loss of contrast in the measured interference 
pattern. 

Thus we have here the important result that the linear part of the potential doesn't 
induce in general any dephasing up to second order in Dyson expansion. In fact we 
show in the next section that dephasing results purely as an effect of the nonlinear 
potential term C 2 Mc 2 A 2 . 



5.2.2. Contribution of the nonlinear part of the potential This calculation requires 
estimating averages of the kind (^4 2 (y, t)A 2 (y' ,t')), which will bring in the second 
order correlation function R" defined in (16). This is straightforward but algebraically 
lengthy. Proceeding in a similar way as done above, exploiting the statistical properties 
(13)-(16) and the already mentioned result I = in relation to (19), then the general 
result for the density matrix and valid up to second order in the Dyson expansion can 
be proved to be: 



,{T) 



'(0) 



32C 2 TT 2 M 2 c 4 A*p xx ,(0) 



dt / dt'R 2 {t-t') 



16tt 2 



dk /dK 



T 

dt 



dt'R(t -t'-k- Ax/c) x R(t — t' — K.- Ax/c) 



(20) 



Remarkably, the second order correlation function doesn't play any role: the first 
order correlation function R(t) and thus the power spectral density S(u>) completely 
determine the system evolution up to second order. Equation (20) implies that the 
diagonal elements of the density matrix are left unchanged by time evolution. This is 
seen by setting Ax = which yields immediately p xx (T) = p xx (0) for every T. 



6. General density matrix evolution for large drift times 



To verify that we have dephasing with an exponential decay of the off diagonal elements 
we need further simplify the result (20) by analyzing its behavior for appropriately 
large evolution times. To this end we start from the following identity 

pT pT -i pT pT poo 

J dt J dt'g{t-t r ) = —jdtjdt'l dLog(uj)e zw{t - t ' ) 
1 f°° 



sin(wT/2) 
ui/2 
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where g(uo) denotes the Fourier transform of the function g(t). Denoting [0, Aw] as a 
frequency interval where g{u>) is slow varying, it is straightforward to show that the 
above identity reduces to 

dt / dt'g{t-t')ttg(0)T, (21) 



for T > (Au;) -1 . Note that for this to happen g(t) doesn't even need being an even 
function. This condition translates what we mean by appropriately long evolution 
time. In section 8.2 we will show that it is equivalent to T > r», where r* is the 
fluctuations correlation time. This is defined below. 

Equation (21) can now be used to evaluate the time integrals appearing in 
(20). This is done by identifying in one case g(t) := R 2 (t) and in the other 
g TT '(t) := R(t + r)R(t + r'), where r and r' stand respectively for — k • Ax and 
— K • Ax, and where the normalized correlation function can be expressed, using the 
generalized W-K theorem (10), as 

R(r) = ^- = 1 /^WS(o,) cos( W r). 

Notice that the integration frequency has a cutoff at tu c — ujp/X, where the Planck 
frequency is ujp := 2n/Tp = 1.166 x 10 44 s _1 . This is consistent with the fact that 
below the scale I = XLp the approximation of randomly fluctuating fields breaks 
down. In alternative this may simply correspond to the fact that the probing particle 
is insensitive to the short wavelengths as a result of its own finite resolution scale -Lr. 
Application of (21) to g(t) := R 2 (t) yields the result 

dt dt'R 2 (t-t') = t*T, (22) 
Jo Jo 

where the correlation time is defined as 

J denoting Fourier transform. On the other hand, application of (21) to g T T'(t) := 
R(t + r)R{t + t 1 ) gives the result 

f dt [ dt'R{t + T)R{t + T r ) =t*T[(J c (t-t')]T, (24) 
Jo Jo 

where the characteristic function T has been defined as 

f" c du uj a S 2 (uj) cosiujt) , s 

r M := J ° re , \' ' ■ (25) 
Jo dww 5 ( w ) 

This is dimcnsionless and satisfies in general the following properties: 

• r(w c t) = r(-cj c t), 

• r(o) = l, 

• r(w c i) < 1, for t ^ 0, 

• r(u} c t) -> 0, for t -» oo. 
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Notice that both the correlation time and the characteristic function T solely depend 
on the fluctuations power spectral density. 

The results (22) and (24) can now be used in equation (20) to yield the neat result 



where 



Pxx' (T) = Pxx' (0) 



F(Ax) := 1 - 



1 - 



32C 2 TT 2 M 2 c 4 Af ) T»T 
h 2 



x F{Ax) 



16tt 2 



dk / dkr[w c (K - k) • Ax/c]. 



(26) 



(27) 



This equation is important and represents one of the main result of this paper. It 
implies that dephasing due to conformal fluctuations dm* indeed occur in general and 
independently of the precise power spectrum characterizing the fluctuations. Without 
the need to evaluate the angular integrals, this follows from the properties of the 
characteristic function T. The fact that r[w c £] < l implies < F(Ax) < 1 with 
(i) F(Ax = 0) = and (ii) F(Ax — > oo) = 1 as special limiting cases. As a 
consequence the diagonal elements are unaffected while the off-diagonal elements decay 
exponentially according to 

'xx'OO - Pxx'(0) 



'(0) 



T 



32C^W^ x ■ 



n 2 



Pxx'(0), 



providing of course that T is small enough so that the change in the density matrix is 
small. Finally, if Sp := p xx >(T) — p xx >(0), we can define the dephasing rate as \8p/po\. 
Thanks to the property F(Ax — > oo) = 1, this converges for large spacial separations 
to the constant maximum value 

32C 2 tt 2 M 2 c 4 A^T 



h? 



(28) 



This result based on the present 3-dimensional analysis of the conformal 
fluctuations can be compared to the analogue 1-dimensional result that PP found 
in [1] . Using a gaussian correlation function from the outset they found 



v^M 2 c 4 A 4 Tg T 



id V2h 2 

where r g stands for some characteristic correlation time of the fluctuations. Identifying 
approximately r* ~ r g , we have (32C27T 2 )/ (-y/ 7r/2) ~ 250, assuming C 2 ~ 1- Thus 
the present 3-dimensional analysis is seen to predict a dephasing rate 2 orders of 
magnitude larger than in the idealized 1-dimensional case. 



6.1. A remark on the validity of the Dyson expansion 
We have found that the change in the density matrix is given by: 

/ Mr 2 \ 2 

In order for the expansion scheme to be effective, the propagation time T must be 
short enough to guarantee that T2[^4o] is small. How short depends of course on the 
statistical properties of the fluctuations, encoded in r„, and on the probing particle 
mass AI. A fullcrcne C70 molecule with (M C70 ~ 10~ 24 kg) gives Mc 2 /h « 10 27 s _1 . 
Therefore the approach is consistent only if correlation time r* , the flight time T and 
field squared amplitude Aq are appropriately small. We will come back on this issue 
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in section 8.2, where it is shown that, in the case of vacuum fluctuations (introduced 
in the next session), it is t* ~ ATp and Aq ~ 1/A 2 . For a flight time T»l ms, typical 
of interferometry experiments, this results in T 2 [^4g] ss 10 7 /A~ 3 . For any reasonable 
value of A > 10 3 the density matrix change is indeed small and the Dyson expansion 
scheme well posed up to second order. In Appendix D we estimate the fourth order 
term in the expansion, which will also yield a term proportional to Aq. It will be 
shown that its contribution in fact vanishes under quite general circumstances. This 
puts the result (26) on an even stronger basis. 

7. Explicit dephasing rate in the case of vacuum fluctuations 

The result (26) is quite general. The only ingredients entering the analysis so far 
have been: (i) a spacetime metric g a b = (1 + A) 2 n a i, with A = 0(e <C 1), (h) a 
randomly fluctuating conformal field A satisfying the wave equation d c d c A = 0, and 
(iii) the isotropic fluctuations characterized by an arbitrary power spectral density 
S(uj). The dephasing then occurs as a result of the nonlinearity in the effective 
potential V = Mc 2 [C\A + C 2 A 2 }. 

A particularly interesting case, potentially related to the possibility of detecting 
experimental signs of quantum gravity, is that in which the fluctuations in A are the 
manifestation, at the appropriate semiclassical scale £ = XLp, of underlying vacuum 
quantum fluctuations close to the Planck scale. Strictly speaking the presence of 
the probing particle perturbs the genuine quantum vacuum state. For this reason it 
would be appropriate to talk of effective vacuum, i.e. up to the presence of the test 
particle. By its nature, the present semiclassical analysis cannot take in account the 
backreaction of the system on the environment. Therefore we simply assume that the 
modifications on the vacuum state can be neglected as long as the probing particle 
mass is not too large and the evolution time short. We thus model the effective 
vacuum properties of the conformal field A at the semiclassical scale on the basis of 
the properties that real vacuum is expected to possess at the same scale. It is a fact 
that vacuum looks the same to all inertial observers far from gravitational fields. In 
particular, its energy density content should be Lorentz invariant. This can obtained 
through an appropriate choice of the power spectrum S(u>). 

7.1. Isotropic power spectrum for 'vacuum' conformal fluctuations 

According to the above discussion we expect the average properties of A above the 
scale £ to be Lorentz invariant. In particular, the interesting quantities derived in 
section Appendix A. 5 



should be invariant. As discussed in Appendix A, for a stationary, isotropic signal, 
the averages (•) can in fact be carried out through suitable spacetime integrations over 
an appropriate averaging scale L 3> I. In alternative they can be expressed as in the 






Dephasing due to conformal fluctuations 



16 



above integrals depending on the power spectrum and adopting a high energy cutoff 
set by k\ := 2it/{\L P ). 

The problem of the Lorentz invariance of the above quantities has been discussed 
in details by Boycr [14] within his random electrodynamics framework. He showed 
that the choice 

S ^ ■= f^< ™ 

guarantees a Lorentz invariant measure d 3 fc/w(fc) (see also [37]) and implies an energy 
spectrum g(u>) oc w 3 , also shown to be the only possible choice for a Lorentz invariant 
energy spectrum of a masslcss field. The combination of the constants h, G and c 
gives the correct dimensions for a power spectrum (i.e. L 3 ), while the factor 1/2 
guarantees that the resulting energy density is equivalent to that resulting from the 
superposition of zero-point contributions fuv/2. A final important point, which should 
not be overlooked, is that Lorentz invariance is preserved provided the cutoff k\ is given 
by the same number for all inertial observers, as also discussed in details by Boycr. In 
other words this means that the critical length that sets the border line between the 
random field approach and the full quantum gravity regime is supposed to be the same 
for any inertial observer. It represents some kind of structural property of spacctimc 
and not an observer dependent property. Accordingly it must not be transformed 
under Lorentz transformations. It is important to note that this requirement will also 
be satisfied when we employ an effective cutoff set by the particle Compton length. 

Using (29) the normalized correlation function can be found explicitly from the 
generalized Wicner-Khintchinc to be 

R(r) := = 2 + ^ff^ \ . (30) 

C LU C T UjfT 1 



Figure 1. Plot of R?(t — t'). The adimcnsional variable cr is basically t — t' in 
units of the correlation time t, . It is seen that r fss r* corresponds to the first of 
the secondary peaks. 
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The peak of the autocorrelation function is linked to the fluctuations squared 
amplitude and gives explicitly: 

C(Q) = A 2 = ^, (31) 

implying (A 2 } = 1/(2A) 2 . The correlation time and characteristic function follow 
from equations (23) and (25) as: 

r* = ^AT P , (32) 

and 

3sin(g) 6cos(cr) 6 sin(q) 

1((J) = h 2 3 ' 

a a" o~° 

where a = uj c t is a dimensionless variable. The plot of the squared normalized 
correlation function R 2 (t — t') is shown in figure 1: t — t' = t* corresponds to the 
first secondary peak in the curve, where the correlation in the fluctuations is reduced 
of ~ 70%. This fully motivates the choice of t* to represent the correlation time. 

The explicit form of the characteristic function can be used in (27) to evaluate 
the remaining angular integrals and find the detailed expression for the density matrix 
evolution valid for all (x,x'). Isotropy implies that the result must depend on |x — x'| 
only. The integration is straightforward and yields the result: 

_. . 3 / sin cr cos cr \ 

F(a):= l-— 2 (l —). (34) 

Substituting the results (31), (32), (33) and (34) into (26) yields the explicit result for 
the dephasing rate, valid for vacuum fluctuations described by S ~ l/u>: 



1 f MY f T\ ^Trlx-x'l , 

3^Up) {^J xF \- L r- L] ' 



Figure 2. Plot of the function F(cr) in the range <r = 0..10, where cr = 
27T |x — x'| / (\Lp). The curve tends very rapidly to the limiting value of 1 and for 
spacial separations |x — x'| which arc slightly larger than i = \Lp the dephasing 
rate converges rapidly to its maximum value. 



0.8 
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where we considered C2 ~ 1 and where 

m p : = ~o7FT = \\-^ = 2-176 x l(T 8 kg = 1.310 x 10 19 amu 
crip V Cr 

is the Planck mass. The function F is plotted in figure 2. It enjoys the properties 
F(0) = and F(a) — » 1 for cr » 1, so that for |x — x'| > 101 the decoherence rate 
converges rapidly to its maximum value. 

8. Discussion 

8.1. Probing particle resolution scale and effective dephasing rate 

Equation (35) is an important result. It gives the dephasing rate in the density matrix 
of a quantum particle propagating in space under the only action of a randomly 
fluctuating potential due to spacetime vacuum conformal fluctuations. The fact that 
it predicts an exponential decay of the off diagonal elements of the system density 
matrix (which is the distinctive feature of genuine quantum decoherence) is interesting 
as a further confirmation that certain effects involving quantum fluctuations can be 
mimicked by means of a semi-classical treatment in the spirit of Boyer [14. 15]. 

A significant feature of our dephasing formula is the quadratic dependence on the 
probing particle mass M, which comes as a consequence of the underlying non linearity. 
The coefficient 1/(3A 3 ) sets the overall strength of the effect. It is proportional to Aq 
and to the fluctuations correlation time r* : the more intense the fluctuations, the larger 
the dephasing and the longer the various directional component stay correlated, the 
higher their ability to induce dephasing. We have found t* w ATp, in such a way that 
the correlation time directly depends on the spacetime intrinsic cutoff parameter A. 
According to this picture all the wavelength down to the cutoff £ = XLp should be 
able to affect the probing particle. However an atom or molecule is likely to possess 
its own resolution scale Lr. Thus, whenever r»c < Lr, the ability of the fluctuations 
to affect the particle would be reduced, as they would effectively average out. To 
characterize this feature of the problem we write, in analogy to £ = XLp, 

Lr = ArLp, 

and use Ar as a new, particle dependent, cutoff parameter. In general it is Ar > A. 
The new effective correlation time is now given by t* f=a ArTp. The distance traveled 
by the fluctuations during a correlation time is L* = cr* = 2Lr/3. Thus the effective 
correlation distance L t basically corresponds to the particle resolution scale: short 
wavelengths that do not keep their correlation up up the scale Lr average out and 
cannot affect the probing particle. The new, effective dephasing rate results by 
substituting A with Ar in (35): 

3 \L R J \Mp) \Tp) \ L R J 

8. 2. Validity of the long drift time regime 

We recall that this results holds for 'long' drift times T, i.e. when T > (Aw) -1 , where 
Aw is an appropriate frequency range over which the Fourier transforms of R 2 (t) and 
R(t + r)R(t + t') vary little. We are now in the position to make this precise and 



Po 
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define clearly the limits of applicability of the theory. To this end we consider the 
Fourier transform of R 2 (oj c t): 

S[R 2 (oj c t)](uj) = —$[R\cj)}{u/u c ), 

with R(o~) given in (30). Its plot is displayed in figure 3. The spectrum falls 
to for uj > 2u) c . The value of the peak at uj = is precisely 47r/3, verifying 
that r* = : 5[i? 2 ((jj c r)](0) = 4tt/(3uj c ). The smaller box shows a zoom of the 
plot in the region a G [0,1/100]: the curve is slow varying in this range since 
S[i? 2 (w c T)](0) = 4tt/3 w 4.19 and $[R 2 (u> c t)] (1/100) « 4.13. Similarly it is possible 
to check that the Fourier transform of R(a + T])R(a + rj), where the adimensional 
parameters rj and rj depend on space direction and locations as r\ := — w c k • Ax/c 
and rj := — oj c K. ■ Ax/c, enjoys a similar property: for every choice of rj and rj' , the 
resulting Fourier transform is slow varying in the range a € [0, 1/100]. Following this 
discussion we chose Aw « [0, w c /100]. We can now quantify the concept of 'long drift 
time' by T > 100/cj c . From uj c = 27t/(ArTp) = 47r/(3r*) this yields the condition 

T > 25 r«. 

8.3. Some numerical estimates and outlook 

In summary we have studied the dephasing on a non-relativistic quantum particle 
induced by a conformally modulated spacetime g a b = (1 + A) 2 n a b, where A is a random 
scalar field satisfying the wave equation. The important case of vacuum fluctuations 
can be characterized by a suitable power spectrum S ~ l/u>. If Lr, = XrLp is 
the probing particle resolution scale, the dephasing rate for |x — x'| 3> £r converges 
rapidly to: 

1 /£p\ 3 f M Y [ZS 



3 \L R J \Mp I \T P 



(36) 
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The effective correlation time of the fluctuations is given by t* ps ArTp. The above 
result holds for 'long' drift times satisfying 

T > (10 - 10 2 )ArT p . 

To conclude we want to give some numerical estimates of the dephasing that 
could be expected in a typical matter wave interferometry experiment, e.g. like 
those described in [9], where fullerene molecules have been employed with drift 
times of the order of T =: T cx « 10~ 3 s. Consider e.g. a Cro molecule with 
M = Mc T0 ~ 1-24 x 10 - 24 kg. In comparison to the Planck units we have: 

T CX « 10 40 T P , M C70 « 10- 17 M P . 

Thus, it is clear that the most critical factor controlling the strength of the effect is 
set by the probing particle mass, together with the effective resolution cutoff scale. 
Using these data in (36) we can estimate: 



10 6 



Ar 



/ 10 6 



VI<WH 

This could be used to estimate Ar if we were able to identify within an experiment a 
residual amount of dephasing that cannot be explained by other standard mechanisms 
(e.g. environmental decoherence, internal degrees of freedom). Figure 4 plots Ar 
against \8p/po\: a dephasing rate due to conformal fluctuations in the range 1% — 0.1% 
would imply a resolution parameter in the range Ar ps 10 3 — 10 4 . This would represent 
a lower bound on Ar , as interferometry experiments will get more and more precise in 
measuring and modeling environmental decoherence. Estimating the present typical 
uncertainty of typical interferometry experiments as \5p/po\ 0.01% we get 

A R > 10 3 , for C 70 . 

We remark that such order of magnitudes estimates are consistent with a small change 
in the density matrix and the second order Dyson expansion approach. 

A value for Ar as small as 10 3 would probably approach the intrinsic spacctimc 
structural limit set by A, i.e. £ = XL p. It is interesting to ask what amount of 
dephasing our model predicts, independently of experimental data. To this end we 



Figure 4. Adimcnsional cutoff parameter as a function of the dephasing rate for 
a C70 molecule with a drift time of 10 — 3 s. 



0.002 0.004 0.006 0.008 0.01 



Dephasing due to conformal fluctuations 



21 



(37) 



need to prescribe theoretically the particle resolution scale Lr. Though no obvious 
choice exists, an interesting possibility would be to set it equal to the particle Compton 
length, i.e. 

T k 

LR ~ Wc 

This choice is obviously Lorentz invariant and also motivated by the fact that 
the Compton length represents a fundamental uncertainty in the position of a 
nonrclativistic quantum particle. Indeed, by Hciscnbcrg uncertainty principle, Ax ss 
h/Mc would imply Ap > Mc, implying an uncertainty in the energy of the same order 
of the rest mass Mc 2 . In such a situation QFT would become relevant. Alternatively 
it can also be argued that wavelengths shorter than h/Mc would have enough energy 
to create a particle of mass M from the vacuum. With this choice, equation (36) 
becomes 

S ± 1 (M\ 5 (T\ 

This can be used to predict the amount of dephasing induced by vacuum conformal 
fluctuations. 

In the case of CVo the Compton wavelength is w 10 _18 m ~ 10 18 Lp, corresponding 
to Ar r* 10 18 . For a propagation time of ~ 1 ms this gives a dephasing rate 

^ (M Cro , lms) «10- 44 , 
Po 

which would be negligible and far beyond the possibility of experimental detection. 
Thus, in order to achieve dephasing rate within the current experimental accuracy, 
much heavier quantum particles are needed. In atomic mass units C70 has a mass 
Mc 7a ~ 10 3 amu. Equation (37) applied to a particle with mass M m 10 11 amu and 
with a drift time T w 100 ms gives the estimate: 

— (10 n amu, 100ms) w 10~ 2 . 
Po 

A drift time of ~ 100 ms could possibly be achieved in a space based experiment. 
On the other hand, the need of a quantum particle as heavy as 10 11 amu poses an 
extraordinary challenge. A possibility would be to employ quantum entangles states. 
This is already being considered in the literature, e.g. in [29], where entangled atomic 
states are studied and suggested as a possible improved probe for future detection of 
spacetime induced dephasing. 

The last important point that needs verification is that the condition (4) gave 
earlier at the beginning of this paper is indeed verified: that was required in order 
for the change in momentum due to the fluctuations to be smaller than the 'main' 
particle momentum p = Mv. It read: e 2 ~ (^ 2 ) "C (v/c) 2 . The field effective mean 
quadratic amplitude interacting with the particle is given by (A 2 ) ~ Ar - . Thus we 
have the condition: 

1 hp v 
Ar L r c 

By using the expression for the Planck length and with Lr given by the particle 
Compton length, this yields a condition on the particle mass M: 

M „ v 
— < 2tt-. 

M P c 
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For typical laboratory velocities v/c ~ 10 -6 and, since Mp w 10 19 amu, this condition 
is met for particle masses up to M « 10 13 amu, including the case of C70 molecules or 
the heavier entangled quantum states discussed above. This limit would be reduced 
for slower particles. 

We conclude by remarking that the theory described in this paper is quite general, 
in the sense that as a starting input it only needs a conformally modulated metric and 
a scalar field satisfying the wave equation. Of course, it is important to identify in 
concrete which theories of gravity can actually yield such a scenario. This important 
problem will be the object of future reports, currently in preparation: in [27] a general 
approach for the study of fluctuating fields close to the Planck scale is introduced and 
the resulting framework applied to standard GR; while [28] will consider more general 
scenarios involving scalar-tensor theories of gravity. 
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Appendix A. Stochastic scalar waves and generalized Wiener-Khintchine 
theorem 

Appendix A.l. General solution to the wave equation 

In this appendix we work out some general results holding for scalar stochastic waves. 
Working in units with c = 1, the solution to the wave equation (— df + V 2 )/(x, t) = 
can be written as 

/(x ' t)= (2^/ d " fc/(k ' t)eik ' X ' 

where 

/(k,i) = J d"*/(x,*)e- ik - x 

and where / <G i 2 (R 3 ). The Fourier coefficients take the general form 

/(k,i) = a(k)e- M + b(k)e tkt 

for some complex functions a(k) and 6(k) and where k := |k|. 
The wave can be written as 

a ( k yfc(k-x-t) + 6 (_ k ) e -ifc(k-x-t)l j 

where dk represent the elementary solid angle in momentum space. Using spherical 
coordinates in momentum space with k(i?, ip, k) and employing the notation a^(k) := 
a(i9, ip; k) and := fe(7r — 19, p + Tt; k), we define the directional wave along k(t9, tp) 

as 

a ic (fc)e 4fcc(li - x/c ^ t) +6_ ii (fc)e- 4 ' £c(ic - x/c - t) j , 



/(x,f) 



(2irY 



dk / dkk 7 



4( k .x/c-i):=— - / dkk 
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where for clarity we have restored the speed of light. The general solution can thus 
be written according to the directional decomposition 

/(x,t) = J dk/ t (k.x/c-t). (A.l) 

Appendix A. 2. Stochastic waves 

In this section we review some elementary properties of stochastic signals of one time 
variable t. 

The fluctuating conformal field at a given space location x represents an example 
of a stochastic process f(t). Its properties can be defined in a statistical sense. If the 
possible values of / obey a probability density Pt(f) the statistical average at time t 
is defined as 

(/(*)) := / Pt(f)fdf. 

The variance is 

<[/(*)] 2 ) := / Pt(f)fdf. 

Denoting the probability density of having the particular outcomes f(t) and f(t') by 
Ptf (/;/'); then the autocorrelation or 2-points correlation function of f(t) is defined 
as 

R(t-t>) := (/(t)/(f)) := J PtAfJ')ffdfdf. 

If the probability of having the value f(t') is completely independent from the previous 
outcome f(t) then ptt'(f,f) = Pt{f)pt'{f) an d the stochastic process described by 
f(t) is said to be perfectly uncorrelated. In this case {fit) fit')) = {f{t)) {f{t')). Higher 
order autocorrelation functions can also be defined. The second order correlation 
function is given by 

R"{t-t') = {[f{t)] 2 [f{t')} 2 ) ■= J PtAf,f')f 2 f' 2 dfdf. 

In the case of conformal fluctuations, we assume that the stochastic process is 
stationary, i.e. all average properties do not depend on time, and isotropic. This 
implies that different directional components have the same statistical properties. 
Finally, by assuming the stochastic process to be ergodic, and since it satisfies the 
wave equation, then statistical averages can be replaced by time or space averages 
taken on any given sample function representing the process. 

Appendix A. 3. Generalized Wiener- Khintchine theorem for a stochastic process 
satisfying the wave equation 

We now derive a generalization of the Wiener-Khintchinc theorem linking the power 
spectrum to the autocorrelation function. We assume that a particular sample function 
representing the stochastic process can be written as in (A.l), in such a way that the 
wave equation is satisfied. Some care must be taken in using the Fourier expansions 
relations of the previous section. Indeed, for a given t, f is not in general a square 
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intcgrablc function belonging to L 2 (R 3 ). To circumvent this problem, given a sample 
function of the process and for any given t, we define 

f L (k,t) := / d".T/(x,i)e- lkx , 

where X>l is a 3-dimensional cubic domain of side L. Then, the function 



f L (x,t) := / d"fc/ L (k,ty k - x , (A.2) 

satisfies the wave equation provided that 

f L (k,t) = a L (k)e~ ikt +b L (k)e ikt . (A.3) 
These expressions will be used later while taking the limit with L — > oo. 



Consider now a complex, ergodic, time-dependent stochastic function /(x, t) in 
an n-dimcnsional space with time t, satisfying the wave equation (<9 2 — V 2 )/(x, t) = 0. 
The autocorrelation function of f{x.,i) for any two events (xi,ti) and (x2,<2) = 
(xi + £, t\ + t) is a function of £ and r given by 

C(tr) = </(x 1 ,t 1 )*/(x 2 ,t 2 )), 

and having the property C(— £, — t) = C(£, r)*. For any fixed choice of Xi and ti it 

also satisfies the wave equation, i.e. ^9 2 — C(£, r) = 0. Assuming that, for any 

t, C(£, r) belongs to L 2 (M 3 ) we have 

C(€, r) = J d n k [a(k)e~ ikT + 0(k)e^] e^. (A.4) 

Notice that a(k)* = a(k), /3(k)* = /3(k), i.e. both a(k) and /3(k) arc real. For r = 
(A.4) reduces to the case of stochastic process of one n-dimcnsional space variable: 

C(€,0) = | d"fc[a(k) +/3(k)] e* k «. 

The standard W-K theorem then implies a power spectrum: 

S(k) = a(k)+/?(k). (A.5) 

To determine a(k) and /3(k) we consider the stochastic process for some fixed 
time to- From equations (A.2) and (A.3) we have 

/ L (x,M = J d n kf L (k,t y k *, (A.6) 

with 

f L {k,t ) = a L (k)e- ikt ° + b L (k)e lkto . (A.7) 

We are thus dealing with one space variable stochastic process for which the usual 
W-K theorem holds. Exploiting the fact that the stochastic process is stationary we 
define mean power spectral density as 

S(k) := gmi /' T dtoi(/ i (k,ior/ i (k ) io)}. 

'■ J 

Using equations (A.6) and (A.7) this reduces to 

Comparing with (A.5) we have 

S^=^(\a^)\ 2 + \b^)\ 2 )- 
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Appendix A. 3.1. Real stochastic process If /(x,i) is real it is easily verified that S(k) 
is even, i.e. 5(k) = S(— k) and that C(£,t) is real. Moreover the power spectrum is 
then given simply by 

In this case we get the generalized Wiener-Khintchine theorem for a real, stationary 
stochastic scalar in the form: 

Cfc, r) = / d "* 5 ( k ) cos(k • | - her), 

where we have restored the speed of light. 



Appendix A. 4- Correlation properties of wave components in different directions 

The results that we have come to establish allow to show that wave components 
traveling in different directions are uncorrelated. Evaluating (/*(xi, ii)/(x2, £2)) using 
equation (A.l) we have 

C(€,r) = (r(x > t)/(x + € > t + r)) = 

J dk y dk' (/?(k-x/c-t)4,([k'-x/c-i] + [k'-£/ C -r])). (A.8) 

Using the Wiener-Khintchine theorem in its form as given by (A. 4), restoring the 
speed of light, swapping k to — k in the second integral and splitting the d n k integral 
in its angular and magnitude parts we can write as well 

C(£,r) = / dkC fc (k.£/c-T), (A.9) 

where we defined the correlation function in the direction k as 
1 r°° 

C £ (k • i/c - r) := — — / dfc k n -\a^ky kc ^l c -^ + /^(fc)^^^], 
v Z7r ) Jo 

with a^(k) := a(k) and f3_^(k) := /3(— k). The two equations (A.8) and (A.9) must 
be equivalent. This implies at once the following equation 

(/?(k • x/c - t) 4,([k' • x/c - 1] + [k' • €/c - r])) = S(k, k') C fe (k • €/c - r) 

or, equivalently, since k • x/c has the dimensions of a time, 

(/?(i)/ fc ,(i + r)} =<5(k,k')C fc (r). (A.10) 

We thus see that the fluctuating field can be resolved into components along different 
directions represented by completely uncorrelated functions of just one time variable. 



Appendix A. 5. Isotropic power spectrum and field averages 

In the relevant case of a real and isotropic signal in 3-dimensional space the spectral 
density must also be isotropic, i.e. 5(k) := S(k). We can then introduce an isotropic 
correlation function as 

C{t):=—— dkk 2 S(k)cos(kcT) (A.ll) 
(2tt) Jo 
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and (A. 10) reads simply (f^(t) fy (t + r)) = <5(k, k') C(r). Using this together with 
(A.l), we have 

</(x,t) 2 ) = ^y dk/ t (x,t) y dk'4,(x,t)^ =c y dk y d^k,^, 

so that 

(/(x,i) 2 >=47rC . 
Using this with equation (A. 11) we have 

1 

(2V) 3 

It is also useful to deduce two expressions for the average of the square of the time and 
space derivatives of the field, as these are directly related to the field energy density. 
Using equation (A. 2) and the ergodic property it is straightforward to show that 



(/(x,i) 2 ) = — 3 / d 3 kS(k). 



(at/)2 )=(2^/ d3fcfc25(k) ' 



(27T) 

and 

;i y/|2 > = (2W d3 " 2s < k >- 

Appendix B. Treatment of the term T4 

We consider the correction term T 4 = — ih (^A — AA^j ip — ihv(A x — AA X ) ip, derived 
in section 2. Writing A, in the isotropic and real case, as: 

Ass /dk / ' \kk 2 a{k)e lk ^*- ct \ (B.l) 



11 



where the upper cutoff is set by the particle resolution scale Lr. The power spectrum 
is basically proportional to the square of the Fourier component a(k). For S(k) ~ 1/k 
we have a(k) ~ k~ 1 ^ 2 , so that the effective coefficient appearing in the expansion is 
k 2 a(k) ~ fc 3 / 2 . Thus the short wavelengths close to the cutoff give the most important 
contribution. For this reason we can approximate the field as 

A^JdkJ^ ' R dkk^ 2 A{k~k A )e lk ^ x - ct \ (B.2) 

where the function A(fc — k A ) is peaked around a typical wave number k A . This can in 
principle be selected in such a way that the average properties of (B.l) are equivalent 
to those of (B.2). Effectively we get: 

dkk A 3/2 e lkA ^ x - ct) , (B.3) 

so that the conformal field is approximated as a fast varying and isotropic random 
signal characterized by a single typical wavelength X A = 2ir/k A . In relation to the 
fluctuations ability to affect the particle, this will be close to the particle resolution 
scale, i.e. we put X A = kLr, with n > 1. From (B.3) we now have: 

A w -ik A cA = cA = —A, 

where we used = h/Mc. The space derivatives yield: 

A tX w ik A 5/2 J dkk x e lkAi ^ x - ct) = 0. 
Using these two relations in (6) yields the result (7). 
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In this appendix we prove that the result 
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I:= [ dt 
Jo 



2 [ dt'f(t -t')- [ dt'f(t - 1' - k • Ax) 
Jo Jo 



= 



holds for an arbitrary even function f and in the limit k ■ Ax/T — > 0. 

For simplicity let A := k ■ Ax. Defining the variable r := t — t' we have 

fdt'f{t-t') = [dTf(r), 
Jo Jo 

while, with r := t — t' — A, 
[ T dt'f(t -t' — A) = f A dr/(r) 

Jo Jt-A-T 

Introducing the primitive of f(t) 



t-A nt-A—T 

dr/(r) - / dr/(r). 
o Jo 



F(t) := / dr/(r) 
Jo 

we can re- write / as 

1 = 2 [ dtF(t)- [ dtF(t-A)+ [ dtF(t- A-T). 
Jo Jo Jo 

Performing a further change of variable z := t — A the second integral reads 
dtF(t-A)= / dzF(z):= / dzF(z)- / dzF(z). 

J-A Jo Jo 

Performing a similar operation on the third integral we obtain 

r-T j-T-A r-A r-A r-T-A 

1 = 2 dtF(t)- dzF(z)+ dzF(z)+ dzF(z)- dzF(z). 

Jo Jo Jo Jo Jo 

Now we use the information T ^ A and f(t) = f(—t). As an elementary consequence 

we have that F(t) is F(t) = —F(—t), implying that 

x /--x 



fX f-x 
/ dTF(r) = / drF(r) 
Jo Jo 



>0 JO 

Approximating T ± A f=a T and swapping the sign of the upper integration bound 
appropriately we have 











[ dtF(t) - 


f dzF(z) - ( dzF(z) 




[ dzF(z)- ( dzF(z) 


'o 


Jo Jo 




Jo Jo 



The integrals from to A can all be neglected exploiting again the fact that T ^> A 
and we obtain 



I « 2 [ dtF(t) -2 ( dtF{t) = 0. 
Jo Jo 

The result in exact in the limit A/T — *• 0. 
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Appendix D. Fourth order term in the Dyson expansion 

The fourth order propagator in the Dyson expansion is 

/ _ ■ \ 4 r r r t^ ,t< 2 > ,t< 3 > 
K 4 (T) := fyj J dt« J dtWJ dt^j dt( 4 )^(i( 1 ))F(t( 2 ))i?(t( 3 ))i?(t( 4 )),(D.l) 

implying a fourth order term in the expression for the density matrix evolution given 

by 

(^-^V^V^'V^'V' J dWxy^l-/ dV 4) |y (4) )(y (4) l x 

x(^(y (1) ,t (1) )ny< 2 \t< 2 ))ny (3) ,i (3) )V(y (4) )i (4) ))p . 

since the potential is V = Mc 2 {C 1 A+C 2 A 2 ) the average yields one term proportional 
d A 4 : 

V(yU , t^)V(yW , ^)) F(y (3) , i (3)) F(y (4) , t (4) } V 

where := <W — yW • kW. This requires knowledge of the 4-point correlation 
function, involving the average of the product of 4 directional components evaluated 
at different points. For a real random process having a zero mean and gaussian 
distribution the 4-points function reduces to [38, 39]: 

(/A^l/A^)) = (/ k (r (1) )/ k ,(r( 2 >)) (4(t( 3 >)4,(t( 4 >)) 

+ (/ k (- (1) )/ k '(- (3) ))(4(- (2) )4'(- (4) )) 

+ (/ k (r (1) )/ k ,(rW))(/ k (r( 2 ))/ k ,(r( 3 ))}. 
We can now use equation (14) to express the 2-point correlations: 
(4(r (1) )/ k ,(r( 2 ))/ k (r( 3 ))4,(r( 4 ))) = 



-5(k( 1 ),k( 3 ))i?(r( 1 )-r( 3 )) ( 5(k( 2 \k( 4 ))i?o 



(2) _ 



This implies that the term T4^4q] deriving from ( K^po) has the structure: 



-I 2 



dk« dk^ / di« / dt^6(k {1 \k^)R(r^ - r< a >) 



UA A 



where all upper bounds in the time integration can be set equal to T by appropriate 
normal ordering [36] . By carrying out one of the two angular integrations we have: 
i- -i ? 

p fT pT 

UAt 

The double time integral can be simplified using the general result (21) and we have 



dkW / dt« / dt&R[tW - t « - k« • (y« - y ^)] 



T 4 [A 4 ] ~ 3 



Jdk (1 ^[R(t + T)} (0)T 
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where ^ denotes Fourier transform and r := — lo 1 ) ■ (y' 1 ' — y^)- The Fourier transform 
can be evaluated using that R(t + r) = C(t + t)/Cq. Then 



5 [fl(t + r)] (c 



dtC(t + T)e" 



da/ / dta/ S(u)')e lu cosa/(i + 7") 



2C (2 

Integrating with respect to t gives 
ff[ii(t + r)] H = 
This vanishes for 
ff[ie(t + r)] (0) = 



C (27rc) 3 y 

Wf Ww ' 25(w,) /I d ' ^ 



i(u) — ui)t iuj' r , i(-~u> —uj) 



t — r 

e J 



2C (27rc) 3 7 
This vanishes for lj > u c . For w = 0we have: 
1 



Auo'uo' 2 S{uj') 8{uo' - u>)e luJ ' T + 8{-uj' - uj)e- w ' T 



2C (2nc)3 J 

Carrying out the frequency integral and using the properties of the 6 function we have: 

1 



$[R(t + T)] (0) 



lim u> S(uj). 



C {2irc) 3 w-o" 

In the interesting case S ~ 1/ui this tends to 0, in such a way that T 4 [A$\ ~> and 
the fourth order Dyson expansion term doesn't give any contribution. 
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